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Abstract 

We study the propagation of Cerenkov photons generated by Very 
High Energy 7— rays and hadrons in the atmosphere. The photon pro- 
duction height distributions are estimated from semi-empirical meth- 
ods and compared with those derived by standard simulation tech- 
niques. Incident spectra at various observation altitudes are then 
derived after applying wavelength dependent corrections due to pho- 
ton attenuation in the atmosphere during the propagation of photons 
from the height of production to the height of observation. These 
are generated both for 7— and hadron primaries of various energies. 
The derived production height distributions agree very well with those 
generated by the simulation package 'CORSIKA' at all energies and 
for both 7— ray and proton primaries. The incident photon spectra 
are found to be both altitude dependent and primary energy depen- 
dent. The peak of the incident spectrum shifts towards the shorter 
wavelength with increasing altitude of observation for a given primary. 
Also the peak of the photon spectrum shifts towards the shorter wave- 
length with increasing energy of the primary at given altitude. The 
fraction of the UV component in the incident Cerenkov spectrum is 
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estimated both for 7— ray and hadronic primaries at various observa- 
tion altitudes and energies. Hadron generated Cerenkov spectra are 
marginally richer in UV light and the difference increases slightly at 
higher altitudes. The fraction of the UV to the visible light in the 
Cerenkov spectrum could be a useful parameter to separate 7— rays 
from cosmic ray background only if one can measure this fraction very 
accurately. 



1 Introduction 

In the highest energy end of the electromagnetic spectrum, the energy range 
between 50 keV and 30 GeV has been successfully covered by satellite based 
missions 0. At higher the energies steeply falling spectrum makes it im- 
possible to observe the 7— ray sky by satellite based detectors. The ground 
based atmospheric Cerenkov technique has proved to be the most sensitive 
technique in exploring the celestial 7— rays in the energy range 300 GeV - 10 
TeV. The basic detection technique in this case is to deploy one or several 
parabolic mirrors at the observation level fitted with fast photo tubes at the 
focus. Very high energy celestial 7— rays initiate an electromagnetic cascade 
in the atmosphere as they enter the atmosphere. The electrons and positrons 
in the cascade, being relativistic, emit Cerenkov light as they propagate down 
the atmosphere resulting in a faint (~ 10~ 4 times the brightness of the star 
light background) flash of light lasting a few nanoseconds. This fast Cerenkov 
flash is detected electronically by coincidence technique. However the main 
drawback of the technique is the presence of much more abundant cosmic ray 
background which severely limit the sensitivity of this technique. Various in- 
genious techniques are employed to distinguish between the Cerenkov light 
produced by cosmic 7— rays from that by the cosmic rays. Among these, the 
imaging technique, wavefront sampling technique, angular and spectral sepa- 
ration techniques are currently employed by different experiments |2|, f|, |||, g]. 

In spite of the necessity of carrying out these observations during clear, 
moon-less nights which scales to only about 10-15% duty cycle, the technique 
has been successful mainly due to two reasons. Firstly, a large collection 
area which is ~ 10 4 — 10 5 times larger than the satellite detectors provides 
a high event rate. Secondly, with high angular resolution and some recently 
developed 7— hadron separation techniques || |J provide high sensitivity 
to the method. Large ground based telescopes based on this technique are 
now being built to cover the intermediate energy range of 30 - 300 GeV 

H n, g7|, gg . 



In the absence of standard mono-energetic beams of cosmic rays or 7— rays, 
one has to depend on the simulation techniques to understand and optimize 
the detector response to these radiations. For this purpose we have carried 
out detailed simulations with VHE 7— rays and cosmic rays of various pri- 
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mary energies initiating extensive air showers in the atmosphere. In order 
to ensure the reliability of the conclusions drawn from these studies it is 
imperative to verify the correctness of the basic results derived from these 
simulations. There are a few simulation packages painstakingly developed 
by high energy physicists and we chose one of them called CORSIKA [|J 
developed by the KASCADE collaboration. 

CORSIKA (version 562), f|, ^] has been used here to simulate Cerenkov 
light emission in the earth's atmosphere by the secondaries of the extensive 
air showers generated by cosmic ray primaries or 7— rays. This program 
simulates interactions of nuclei, hadrons, muons, electrons and photons as 
well as decays of unstable secondaries in the atmosphere. It uses EGS4 
code || for the electromagnetic component of the air shower simulation and 
GHEISHA[|K]] for the simulation of hadronic interactions at TeV energies. 
The Cerenkov radiation produced within the specified bandwidth (300-550 
nm) by the charged secondaries is propagated to the ground. The US stan- 
dard atmosphere parameterized by Linsley [|I^] has been used. The posi- 
tion, angle, time (with respect to the first interaction) and production height 
of each Cerenkov photon hitting the detector at the observation level are 
recorded. 

We have mainly used Pachmarhi (longitude: 78° 26' E, latitude: 22° 
28'N and altitude: 1075 m) as the observer's location where an array of 
25 Cerenkov detectors each of area 4.35 m 2 is deployed in the form of a 
rectangular array. We have assumed 17 detectors in the E-W direction with 
a separation of 25 m and 21 detectors in the N-S direction with a separation 
of 20 m. This configuration, similar to the Pachmarhi Array of Cerenkov 
Telescopes (PACT) but much larger, is chosen so that one can study 
the core distance dependence of various observable parameters. Simulations 
were also carried out at sea level and a location at an altitude of 2 km above 
mean sea level. Mono-energetic primaries consisting of 7— rays and protons 
incident vertically on the top of the atmosphere with their cores at the centre 
of the array have been simulated in the present study. 

The longitudinal shower development profiles of Cerenkov photons in the 
atmosphere have been generated using CORSIKA for 7— ray and proton pri- 
maries of various energies. Tracks of charged particles are followed through 
in steps 20 g cm~ 2 [III, EH]. 



In the present work the Cerenkov photon growth curves have been ob- 
tained from the simulations for proton and 7— ray primaries of various en- 
ergies and for the three observation levels mentioned above. These have 
been compared with those derived from independent empirical relations. The 
Cerenkov photon spectra as seen at the observation levels are then derived 
for various primary energies of 7— rays and protons after taking into account 
of the wavelength dependent attenuation in the atmosphere. Then we esti- 
mate the relative fraction of UV photons in these spectra as a function of 
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primary energy and observation level. 



2 Calculation of longitudinal profiles 

2.1 7— ray primaries 

For showers initiated by photons of energy E 1 the longitudinal shower de- 
velopment curve is derived from the equation for cascade curve under ap- 
proximation A. The average number of electrons of all energies N(E 1 , x) as a 
function of atmospheric depth x, is given approximately, in the region where 
the number of particles is large, by E 
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7V(£: 7 ,x) = == exp (t(l - 1.51ns)) 



>»(t) 



where t is the depth x measured in units of radiation length in the atmosphere 
(37.2 g cm~ 2 ) and E t is the the electron critical energy at the depth t. This is 
the energy below which the electron multiplication stops. It is approximately 
the energy at which radiation losses and collisional losses are equal and has 
a value 84.2 MeV in air. The shower age parameter, s, which is a measure 
of the stage of the longitudinal shower development, is given by: 

3t 



(* + 21nf) 

s = 1 at shower maximum while s < 1 above the shower maximum and and 
s > 1 below. The relation between the atmospheric height h, (as measured 
from sea level in m) and the depth, x (g cm~ 2 as measured from the top of 
atmosphere) is taken from Rao| 



m 



h = (6740 + 2.5a;) In C-^~) 



where the scale height is dependent on the atmospheric depth. The refractive 
index n at a given height h is assumed to be: 



n = 1 + n 

where n is given by: 

n = 2.9 x 10~ 4 e~7T5o 

For calculating the electron growth curve the atmosphere is divided into 
slabs of thickness 333 m, after the height of first interaction of the primary. 
The electron threshold energy for the production of Cerenkov radiation at a 
depth t is given by 
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E t h = i - 1 m e 

^-(l+no)- 2 J 

where m e is the electron rest mass (0.51 MeV). The fraction of electrons with 
energy above the Cerenkov threshold is given by [|l(J: 

/0.89£ - l.2\ s , . , 

fc = ( -p ° (1 + 0.0001^)- 2 

where £ = 26 MeV when s < 0.4 

= 44 - 17(s - 1.46) 2 Mel/ for s > 0.4 

The number of electrons above the Cerenkov threshold is then computed 
in each of the slabs. For each straight section of the electron path, the number 
of Cerenkov photons produced per unit path-length is given by, 

dN ( 1 1 \ ( 1 
6.28a — 



dl ' VAi X 2 J V ft 2 ™ 2 ) 

in the wavelength band bounded by Ai and A 2 and ft is the electron velocity 
and a is the fine structure constant. Assuming all the electrons in the cascade 
to be close to the shower axis and they travel with the velocity close to that of 
light in vacuum, the number of Cerenkov photons produced in a path-length 
dl (= 333 m in the present case) and in the wavelength range 300-550 nm at 
a depth x g cm~ 2 is then given by, 

Nph = 0A02N ( dl 
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2.2 Proton primaries 



Because of the differences in the kinematics of shower development in the case 
of hadronic primaries the longitudinal development profiles too are different 
as compared to that of 7— ray primaries. Using the scaling model for nuclear 
interactions, Gaisser & Hillas |T7| find that the average number of particles, 



N(E P , x) at an atmospheric depth of x g cm~ 2 , measured from the first point 
of interaction, in a shower initiated by a proton of energy E p (GeV) can be 
adequately represented by the empirical relationship |?D], O] 

E / t \ tm 
N(E p ,x) = S ^e t -(Pj e-t 

where 

t m (E p ) = jln 

where 5*0=0.045, e=0.074, xq= 37.2 g cm~ 2 , t — f and A =proton interaction 
mean free path in air (70 g cm~ 2 ). We chose the electron energy spectrum 
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Figure 1: Cerenkov photon growth curves in the atmosphere for 7— ray of 
(a) 500 GeV & (b) 1 TeV and protons of (c) 1 TeV & (d) 2 TeV primaries. 
The histograms indicates the simulation results from CORSIKA while the 
smooth curves indicates the results from analytical calculations. The agree- 
ment between the two curves is good. 

given by Zatsepin and Chudakov which is independent of the stage of cascade 
development S HI- 

7 c 

F{E,x) dE = 0.75 N(x)^-^ 

where F(E,x) and N(x) are the number of electrons of energy E and all 
energies respectively at a depth x, 5 = 2.3E/k and k = 72 MeV with E in 
MeV. Even though this approximates well at shower maximum in the energy 
range 20-300 MeV, the minor departure at other depths affects only the 
shape of the lateral distribution of Cerenkov photons and does not affect the 
present estimates significantly. The integral electron energy spectrum is then 
derived: 

F(E,x) = 0.75 N(x) In (l + ^—E^j 

Substituting E t h for E, the number of electrons above the Cerenkov 
threshold at a particular depth are calculated using the above equation. 
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Table 1: Quantitative comparison of the simulated and calculated Cerenkov 
photon growth curves in the atmosphere for 7— ray and proton primaries at 
two primary energies. 



Results 


type of 


Energy of 


m j 1 it r 

Total # of 


from 


primary 


primary 


Cerenkov 






(GeV) 


photons (xlO 7 ) 


Simulations 


7— rays 


500 


2.25 


Protons 


1000 


2.35 


7— rays 


1000 


4.55 


Protons 


2000 


4.81 


Calculations 


7— rays 


500 


2.22 


protons 


1000 


2.35 


7— rays 


1000 


4.56 


protons 


2000 


4.72 



The longitudinal Cerenkov photon profiles are then derived for 7— rays 
of primary energy 50, 100, 250, 500 and 1000 GeV and protons of energy 
100, 200, 500, 1000 and 2000 GeV. The growth curves are simulated using 
CORSIKA for two primary energies for 7— rays (500 GeV & 1000 GeV) and 
protons (1000 & 2000 GeV). Figure shows a comparison of the simulated 
longitudinal profiles for 7— rays (a & c) as well as protons (b & d) with 
calculated profiles. The simulated profiles (histograms) are averaged over 30 
showers for 7— rays and 50 showers for protons. The agreement between the 
simulated and calculated longitudinal shower development profiles in terms 
of Cerenkov photons is good for 7— ray as well as proton primaries of energies 
considered here. Table |l] summarizes the total number of Cerenkov photons 
in these showers. The total number of Cerenkov photons agree with that from 
simulations within 1 & 2% for 7— ray and proton primaries respectively. 

The depth of shower maximum is expected to increase logarithmically 
with primary energy both for electromagnetic as well as hadronic showers 
as more cascade generations are required to degrade the secondary particle 



energies p0| , |TTf . Figure |2| shows a plot of the depth of shower maximum as 
a function of primary energy, both for 7— ray and proton primaries in order 
to illustrate this point as well as to confirm the correctness of the present 
calculations. The energy dependence of the depth of shower maximum is 
parameterized as follows: 

Xmax = 30.3 ln-E-y + 118.8 g arT 2 

for 7— rays and, 

Xmax = 32.6 In E p + 118.6 g cm~ 2 

for protons. 
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Figure 2: The variation of the depth of shower maximum with primary 
energy for 7— rays and protons. 
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3 Photon attenuation in the atmosphere 



To study the attenuation of Cerenkov light in the atmosphere Elterman's 
atmospheric attenuation model p2 is used, which provides the attenuation 



coefficients for the Rayleigh and aerosol scattering as well as ozone absorption 
in an altitude dependent form for the wavelength range 270-1260 nm [^2[ . 

Atmospheric attenuation model computes optical parameters spectrally 
and with altitude as follows: (1) pure air attenuation parameters are deter- 
mined by utilizing Rayleigh scattering cross sections with molecular number 
densities from standard atmosphere; (2) ozone absorption parameters are 
derived from coefficients applied to a representative atmospheric ozone dis- 
tribution; (3) seven sets of aerosol measurements are compared and a profile 
of aerosol attenuation coefficients as a function of altitude is developed. Tab- 
ulation permits calculations for vertical path transmission at one kilometer 
intervals up to an altitude of 50 km, individually for each attenuating com- 
ponent or for overall atmospheric extinction (molecular + ozone + aerosol). 
Above 700 nm light of night sky (LONS) increases rapidly because of emis- 
sion lines of OH and water bands in upper atmosphere while the intensity 
of Cerenkov light drops considerably( inversely proportional to the square 
of wavelength). At lower wavelengths (< 300 nm) light undergoes strong 
absorption by ozone molecules in the upper atmosphere. 

The number of photons, integrated over the bandwidth of 300-550 nm 
(dictated by the photo-tube band-width) and transmitted through a slab (1 
km thick) using the corresponding extinction coefficient is given by: 

N = N ie - B 

where N is the number of photons transmitted through the slab with ab- 
sorption coefficient B and iVi is the mean number of photons entering the 
slab. For each altitude, the photon transmission spectrum is convolved with 
the Cerenkov photon emission spectrum to get the transmitted photon spec- 
trum. Showers pass through several such slabs with different transmission 
coefficients, resulting in a modified longitudinal development profile for the 
Cerenkov photons that reach the observation level. The ratio of the total 
number of photons within the band-width in a shower received at a particular 
observation level to the total photons produced is defined as the transmission 
coefficient. The average transmission coefficients (T c ) are listed in tables ||] 
and [| for 7— ray and proton primaries respectively for Pachmarhi altitude. 
T c can be expressed as a power law in primary energy, E as 

T c = aE b (1) 

The values of a and b are listed in table f| both for 7— ray and proton 
primaries at three different observation levels. The variation in the average 
transmission coefficients for Cerenkov photons in the atmosphere is shown 
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Figure 3: Average transmission coefficient for Cerenkov photons in the at- 
mosphere for primary (a) 7— rays and (b) protons of various energies. The 
straight lines show the fits. The three plots in each panel correspond to three 
observation levels (1) sea level, (2) Pachmarhi (1 km a.s.l.) and (3) 2 km 
(a.s.l.) altitude. 

in figure for (a) 7— ray and (b) protons as a function of primary energy at 
three different observation levels, illustrating the power law behavior. 

The increase in T c with primary energy implies that higher energy pri- 
maries penetrate deeper in the atmosphere and hence pass through lesser 
air mass. The proton primaries reach the shower maximum lower down in 
the atmosphere compared to a 7— ray primary of the same energy since the 
interaction mean free path in air for the former is nearly twice the radiation 
length. As a result the protons have marginally larger (~ 1.2% at 100 GeV 
to ~ 0.6% at 1 TeV) average transmission coefficient compared to a 7— ray 
primary of the same energy. 

The attenuation of optical photons due to Rayleigh and aerosol scattering 
is more significant at lower altitudes. As a result the Cerenkov photon trans- 
mission is expected to depend on the observation level as also seen in figure 
H The average photon transmission coefficient increases almost linearly with 
decreasing atmospheric depth and the rate is ~ 0.1% for every 100 g cm~ 2 
of air. 

Figure f| shows the longitudinal development profiles of Cerenkov photons 
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Table 2: Shower size and average transmission coefficients corresponding 
to Pachmarhi altitude for primary 7— rays of various energies (band-width: 
300-550 nm). 



Primary 


#ofC 


Average 


Depth of 


energy 


Photons 


Transmission 


shower max. 


(GeV) 


produced 


transmitted 


coefficient 


(g cm~ 2 ) 


50 


2.02 x 10 6 


1.29 x 10 6 


0.641 


241.3 


100 


4.16 x 10 6 


2.70 x 10 6 


0.649 


254.0 


250 


1.08 x 10 Y 


7.14 x 10° 


0.660 


281.4 


500 


2.22 x 10 7 


1.48 x 10 7 


0.669 


311.5 


1000 


4.54 x 10 7 


3.07 x 10 7 


0.677 


327.6 



Table 3: Shower size and average transmission coefficients corresponding 
to Pachmarhi altitude for primary protons of various energies (band-width: 
300-550 nm). 



Primary 


#ofC 


Average 


Depth of 


energy 


photons 


transmission 


shower max. 


(GeV) 


produced 


transmitted 


coefficient 


(g cm~ 2 ) 


100 


2.17 x 10° 


1.42 x 10 6 


0.657 


267.4 


250 


5.56 x 10 6 


3.70 x 10 6 


0.665 


296.1 


500 


1.15 x 10 7 


7.74 x 10 6 


0.673 


327.6 


1000 


2.37 x 10 7 


1.61 x 10 7 


0.681 


344.5 


2000 


4.86 x 10 7 


3.34 x 10 7 


0.688 


362.1 



Table 4: Fitted coefficients a and b in equation 1 for 7— ray and proton 
primaries at three different observation levels. 



Observation 


Fitted 


7— rays 


Protons 


level km a.s.l. 


coefficients 









a 


0.475 


0.487 


(Sea level) 


b 


0.018 


0.019 


1 


a 


0.598 


0.608 




b 


0.018 


0.017 


2 


a 


0.646 


0.697 




b 


0.018 


0.012 
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Figure 4: A comparison of the longitudinal development profiles of Cerenkov 
photons at Pachmarhi level with and without atmospheric attenuation cor- 
rection, a & b are for 7— ray primaries of energy 50 GeV and 1 TeV respec- 
tively while c & d are for proton primaries of energy 100 GeV and 2 TeV 
respectively (band- width: 300-550 nm). 
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Figure 5: Radial dependence of the Cerenkov photon transmission coefficients 
in the atmosphere for (a) 7— ray and (c) proton primaries of two different 
primary energies as shown. Also shown are the radial dependence of average 
production heights for each of the primaries (b & d respectively) and their 
energies mentioned above. Simulation results are averages over 43 & 17 
7— ray showers of primary energy 500 GeV and 1 TeV respectively while 47 
showers are used for 1 & 2 TeV proton primaries. 



at production in the atmosphere and that for the photons detected at the 
observation level (i.e. "without" and "with" atmospheric attenuation cor- 
rection) for the same wavelength band. The differences in the two profiles 
is due to the wavelength dependent absorption of Cerenkov photons in the 
atmosphere. From the plots it is clear that the shape of the longitudinal 
profile remains largely unchanged (except for the total number of photons). 
The range of atmospheric heights from which the detected photons originate 
also remains unchanged as well, irrespective of the primary energy or species. 
Thus the height of shower maximum remains unchanged due to atmospheric 
attenuation. 

The average path-length of Cerenkov photons reaching different core dis- 
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tances differ resulting in different attenuation. Figure ^ shows such a radial 
variation of the transmission coefficient for Cerenkov photons generated by 
7— rays and protons of two different energies. Also shown are the radial vari- 
ation of mean production heights. It can be seen that the radial variation 
of transmission coefficients is rather small but varies inversely as the mean 
production height. 



4 Cerenkov photon spectrum 

The Cerenkov photon spectrum at the observation level is an important in- 
put for the design of an atmospheric Cerenkov experiment. The fraction of 
the photon spectrum bracketed by the photo-multiplier bandwidth has an 
important bearing on the sensitivity of a TeV 7— ray telescope. Hence we 
computed the photon spectrum at different observation levels for 7— ray and 
proton primaries of various energies. The bandwidth considered in our cal- 
culation is 270-550 nm. Figure |E| shows photon spectra, corresponding to 
Pachmarhi level, generated by (a) 50 GeV 7— rays & 100 GeV protons (* & 
diamond respectively) as well as (b) 1 TeV 7— rays & 2 TeV protons. It can 
be seen from the figure that the wavelength at peak intensity is a function 
of the primary energy in both the cases. It shifts from around 350 nm at 
50 GeV to around 330 nm at 1 TeV. This is because the higher energy pri- 
maries reach the shower maximum lower down in the atmosphere compared 
to lower energy primaries. As a result, absorption of shorter wavelength pho- 
tons which is primarily due to atmospheric Ozone is comparatively larger for 
lower energy primaries. 



4.1 Fraction of UV component at various altitudes 

The UV filters were used in atmospheric Cerenkov experiments [^| whose 
purpose was two-fold: firstly, it helps reduce night sky background while 
minimizing the loss of Cerenkov light in the blue and near UV range. This 
would provide better stability at higher photo-tube gains which in turn will 
improve the sensitivity of an atmospheric Cerenkov telescope. Secondly, it 
could serve to identify showers with a larger UV fraction in the Cerenkov 
light. Since the Cerenkov light generated by proton primaries traverse lesser 
air mass as compared to 7— ray primaries it is expected to have a larger UV 
content. This property could be exploited to discriminate against hadronic 
showers [p4 |. 

Hence the UV fraction in Cerenkov light at observation level is a useful 
signature. So we estimated this fraction for primaries of various energies. We 
divided the Cerenkov photon spectrum as seen at the observation level into 
two groups viz., the UV range comprising wavelengths 270-300 nm and the 
visible range comprising the wavelength band 300-550 nm. The ratio of the 
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Figure 6: Typical Cerenkov photon spectra at Pachmarhi level, (a) shows the 
spectra from 7— ray primaries of energy 50 GeV and 100 GeV protons while 
(b) shows those for 1 TeV 7— ray & 2 TeV proton primaries. The symbols 
used are * for 7— rays and diamond for protons. 
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Table 5: Ratio of UV (270-300 nm) to visible (300-550 nm) component in 
the Cerenkov light detected at three observation altitudes as a function of 
7— ray and proton energies. 



Primary 


Sea 


level 


1 km 


2 km 


energy (GeV) 


R-f 


R p 


R~/ 


R P 


R-f 


R P 


E 7 /Ep 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


50 /100 


8.23 


9.57 


8.27 


9.67 


8.29 


9.77 


100/200 


8.26 


9.59 


8.29 


9.69 


8.33 


9.80 


250/500 


8.29 


9.60 


8.32 


9.71 


8.37 


9.85 


500/1000 


8.31 


9.61 


8.34 


9.72 


8.41 


9.88 


1000/2000 


8.32 


9.61 


8.36 


9.73 


8.44 


9.91 



number of photons in the UV range to that in the visible range is defined 
as R p for proton primaries and R 1 for 7— ray primaries. Table ||] summarizes 
the R 1 and R p values for different primary energies as seen at three different 
observation altitudes. 

Figure ^ shows a plot of R p and i? 7 as a function of primary energy 
for three different observation levels. Also shown in figure [7] are the relative 
excess of UV content in proton showers vis-a-vis 7— ray showers, as a function 
of primary energy at three different altitudes. 

It can be readily seen from the figure that sensitivity to UV excess is 
better at higher altitudes as compared to sea level. 

Figure |8] shows radial variation of the ratio of UV to visible photons for 
7— ray and proton primaries at an altitude of 1 km. The relative UV fraction 
is marginally higher around the hump region since the photons in this region 
are contributed by higher energy electrons which have shortest path-length 
in air. The radial dependence of the UV fraction is more pronounced for 
proton primaries for reasons already mentioned before. 



5 Discussion and Conclusions 

In spite of the complexity of the interaction kinematics of high energy cos- 
mic rays in the atmosphere, the average behavior as derived from detailed 
simulation studies agrees reasonably well with analytical calculations. This 
demonstrates that the simulation package does take into account almost all 
the interaction characteristics giving credence to the conclusions drawn from 
simulation studies. The position of shower maximum obtained by us agrees 
well with those of Miller & Westerhoff f2§ . 



Armed with the above result we proceed to make analytical calculations 
for the production of Cerenkov light produced by both 7— rays and protons 
of various energies at three different observation altitudes. This is ideal to 
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Figure 7: Variation of the ratio of UV to visible fraction in Cerenkov light 
generated by protons (*) and 7— rays (+) as a function of primary energy, 
at three observation levels: (a) sea level (b) Pachmarhi and (c) 2 km a.s.l. 
(d) The relative excess in the UV content in hadronic primaries as a function 
of primary energy expected at three different observation levels, in the same 
order from bottom to top. 



17 



0.115 



0.105 



m 
c 




r 

a 

« 0.095 

'(/) 
> 

> 
D 



0.085 



0.075 



1 500 GeV 7-rays 

2 1 TeV y-rays 
TeV Protons 

4 2 TeV Proions 




100 200 
Core distancefm) 



300 



Figure 8: Radial variation of the ratio of UV to visible photons in Cerenkov 
light generated by protons (1 & 2) and 7— rays (3 & 4) at two primary 
energies. Alternate error bars are shown for clarity. 
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study the average shower properties as this is much faster than the detailed 
simulations. We then apply a detailed wavelength dependent corrections 
due to the attenuation of Cerenkov photons during their propagation in the 
atmosphere. This has been done for 7— ray and proton primaries of various 
energies and for three observation levels. It has been found that the fraction 
of the transmitted Cerenkov photons at an observation level of 2 km above sea 
level is about 36% more than that at sea level. This is primarily because the 
average distance to shower maximum is reduced resulting in lesser absorption. 
Even though the effective collection area is reduced at a higher altitude, the 
light intensities are larger near the shower core. In addition, increased UV 
content in the Cerenkov light makes higher observation levels better suited 
for atmospheric Cerenkov work. Simulations carried out for Mt. Hopkins 
observatory site (altitude of 2.3 km) suggest an advantage factor in the range 
1.2-1.4 compared to sea level [p9fl , is consistent with the present result. We 
also derived the observed Cerenkov photon spectrum for a 7— ray of energy 
3 TeV at an observation altitude of 2 km. This spectrum agrees well with 



that reported by Mirzoyan, et a/.[|3D| for HEGRA observation level. 

The use of a UV filter to preferentially accept Cerenkov light has been 
suggested even though its efficiency in improving the signal to noise ratio 



of an atmospheric Cerenkov telescope is in doubt ||29[| . On the other hand 



pushing the sensitivity to near UV wavelengths does improve the system 



sensitivity P7[ . Also using suitable filters to suppress the photo-tube 
sensitivity in the visual band, existing imaging telescopes have been modi- 
fied so that observations could be made during moderately moonlit nights. 
This technique has been shown to be successful in detecting TeV 7— ray sig- 
nal from Crab and Mkn 421 despite higher energy threshold and reduced 
sensitivity |25|, 



It has been observed by Zyskin, al, |24| that the relative signal strengths 
in the UV range (200-300 nm) to that in the visible range (300-600 nm) 
increases with the angular size of the image in the visible. It increases up to 
(8 ± 1.8)% for large images {-s/aJ) ~ 0.5° where a and b are the semi-major 
and semi-minor axis of the image in the optical). This could be due to the 
fact that those showers whose maximum is closer to the observation level 
produce larger Cerenkov images. This is equivalent to raising observation 
level which effectively reduces the attenuation of UV photons as discussed in 
§4.1 From the present calculations the relative signal strength at a primary 
proton energy of 4 TeV (energy threshold of Zyskin et al., |24T| ) is estimated 
to be ~ 10% which is quite consistent with their measurement. It may 
be noticed that the wavelength ranges used by Zyskin et al., is broader 
in UV and while the visible band in our estimate was modified to match 
theirs. However the observed Cerenkov photon spectrum falls steeply below 
~ 300 ram, as could be seen in figure |6|. Extending the wavelength range 
below 270 nm does not change the UV content significantly. 
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We find that the relative strengths of UV photons to the visible is higher 
in the case of proton primaries by about 16% at 50 GeV and decreases to 
12% at 1 TeV. Hence, the hadron discrimination efficiency based on the UV 
content in Cerenkov light is relatively better at lower primary energies and at 
near core distances. Thus, measurement of relative UV content of a shower 
could be a good parameter in order to discriminate against hadrons especially 
for large ground based arrays with low energy thresholds (~ 20 — 50 GeV) 
and employing the wavefront sampling technique |35|, |36], [38], ^] . However 



one has to keep in mind that for this technique to be successful one has to 
make a very accurate (better than ~ 1%) estimate of the UV and visible 
light contents of the shower. 
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